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Abstract 

Drop coalescence is central to diverse processes involving dispersions of drops in 
industrial, engineering and scientific realms. During coalescence, two drops first touch 
and then merge as the liquid neck connecting them grows from initially microscopic 
scales to a size comparable to the drop diameters. The curvature of the interface is 
infinite at the point where the drops first make contact, and the flows that ensue as the 
two drops coalesce are intimately coupled to this singularity in the dynamics. Conven- 
tionally, this process has been thought to have just two dynamical regimes: a viscous 
and an inertial regime with a crossover region between them. We use experiments 
and simulations to reveal that a third regime, one that describes the initial dynamics 
of coalescence for all drop viscosities, has been missed. An argument based on force 
balance allows the construction of a new coalescence phase diagram. 

The collision and coalescence of water drops, so essential to raindrop growth and the 
development of thunderstorms, have captivated the attention of the atmospheric science 
community since the early studies by Benjamin Franklin and Lord Rayleigh [28]. Coa- 
lescence also plays a central role in industrial processes involving emulsions or dispersions 
[171 127] . For example, in the petroleum industry, coalescence occurs during dispersed water 
removal and during oil desalting [16] . It is a dominant process in determining the shelf life of 
emulsion-based products such as salad dressing and mayonnaise [21], and it occurs in dense 
spray systems and combustion [3] . Also, sintering of two spherical particles closely resembles 
the coalescence of two dispersion drops in an emulsion [12]. Moreover, the controlled coales- 
cence of drops in microfluidic devices promises a host of potential applications in chemistry, 
biochemistry, and materials science [2]. 

The initial dynamics of coalescence are expected to be universal. The expansion of the 
liquid neck connecting two drops is controlled by the Laplace pressure, which diverges when 
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Figure 1: Coalescence of silicone oil drops with viscosities \i = 58,000 cP (Oh = 440) and 49 cP (Oh = 
0.32). (A) Two pendant drops at the instant they contact, t = 0. (The two bright spots are from back- 
lighting.) We subtract an image taken after the neck has grown to a size r m i n = 0.25^4 from the one at r = 
for (B), fj, = 58,000 cP and (C), 49 cP. 

the curvature of the liquid interface is infinite at the point where the drops first touch. Thus 
the change in topology, as two drops become one, is inextricably linked to a singularity in 
the dynamics. Different regimes of coalescence have been studied [191 EQl I2H EIHl 12S1 DSl [311 
1321 1301 [TJ E21 HI [13 El El ESI E3] . The understanding that has emerged is that coalescence has 
just two dynamical regimes with a crossover region between them: a viscous regime, which 
always dominates at sufficiently early times when the neck radius is microscopically small, 
and an inertial regime that occurs at late times if viscous effects become negligible. We use 
experiments and simulations to show that a third regime, one that describes the true initial 
dynamics of coalescence, has been missed. We present our results in terms of a new phase 
diagram of coalescence that shows the three distinct dynamical regimes. 

Results and Discussion 

In the experiment, two pendant drops with radii A « 0.1 cm are suspended as in Fig. 
[TJA. from nozzles and slowly translated until they touch at their equators. Except where 
otherwise stated, the drops are silicone oil (surface tension 7 = 20 dyn/cm and density p = 
0.97 g/cm 3 ). The liquid viscosity, /1, or equivalently the dimensionless Ohnesorge number, 
Oh = fi/ y/pjA, is varied. We use a high-speed camera and electrical resistance measurements 
[26] to capture the dynamics. 

In the simulations, two isolated spherical drops of radius A are connected by a small neck 
of radius r min = 0.001 A. The dynamics that ensue are determined by solving the full Navier- 
Stokes equations by a finite-element algorithm that we have previously used successfully to 
study diverse situations involving drop breakup [71 [29]. Creeping- flow simulations, where 
inertia is neglected, are also performed. 

Fig. [TJA. shows two drops at the instant of contact, r = 0. To distinguish local versus 
global motion during merging, we subtract an image taken after the neck has grown to a 
small size from the one at r = as shown in Fig. [T]B and C for fluid viscosities \i = 58,000 
cP (Oh = 440) and 49 cP (Oh = 0.32) respectively. At high viscosity, the two drops move 
together rigidly whereas at lower viscosity, the only appreciable motion occurs near the 
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widening neck. This qualitative difference heralds the existence of a new regime. 

The transition between these regimes can be understood from a force-balance argument 
that is based on the observation that in the perfectly viscous (i.e., Stokes) regime, the 
drops outside the immediate vicinity of the neck are rigidly translated towards each other. 
This was shown by Hopper [T9l [20] in an exact analytic solution of coalescence in two 
dimensions (2D); this global motion was also seen in 3D Stokes simulations [21] and the 
early-time asymptotic behavior was later analytically extended to three dimensions (3D) 
[T5] . To be in the Stokes regime, therefore, the force of the neck pulling the two drops 
together must be sufficiently large to produce the required center-of-mass acceleration of 
each drop (i.e., the fluid on each side of the z = plane). The asymptotic acceleration [19] 
is a c . .m. = -y 2 [\n(r min /8A)] 2 /2ir 2 fi 2 A. 

The coalescence is driven by surface tension. Therefore, an upper bound for the inward 
force of the neck on the drops is given by the surface-tension force, F y , around the circumfer- 
ence of the neck at its minimum radius. In 3D: F 7 = 2ir / ~fr m i n . If F y is too small to translate 
the drops, each having a mass m = ^iiA 3 p, the flows cannot be in the Stokes regime. There- 
fore, the Stokes regime can only be achieved when F 7 > ma c . . m . leading to the threshold 
criterion for entering the Stokes regime: 
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Therefore, for 3D drops of finite viscosity, the asymptotic dynamics, in the limit r m i n /A — > 0, 
can never be in the Stokes regime. Below this threshold, F y is balanced by the inertia of the 
drops, ma c o m and the dynamics are governed by local flows. 

An analogous argument in 2D (where F 7 = 27 independent of r m , n , and m = 7iA 2 p) 
suggests a phase-boundary for 2D drops: Oh oc | ln(r min /8A)|/v47T- Thus the logarithmic 
divergence of a c , , m , also precludes the Stokes regime in 2D until r min /A grows to a sufficient 
size. 

Fig. |2jA. presents a phase diagram of the coalescence regimes for 3D drops, which includes 
the inertial and Stokes regimes in addition to this new "inertially-limited-viscous" regime. 
In this regime, inertia and viscosity play a role in the dynamics; the inertia associated 
with each drop moving as a rigid object precludes the system from being in the Stokes 
regime. The inertially-limited-viscous to inertial crossover was previously determined [26] to 
be r min /A oc Oh, in contrast to earlier work that had suggested that this crossover occurs 
at r m inj A oc Oh 2 [T5], |30j [T]. (Previously, there were believed to be only two coalescence 
regimes — a viscous one and an inertial one — so this crossover is referred to as the viscous- 
to- inertial crossover in the literature.) The inertially-limited-viscous to Stokes transition is 
described by eqn. [TJ We emphasize that, contrary to earlier studies, we find the Stokes and 
inertial regimes do not share a phase boundary; they are both preceded by the inertially- 
limited-viscous regime. Thus, at early times, a model of pure Stokes flow for the coalescence 
of spheres is never valid. We note that this is reminiscent of the singularity in drop break-up 
[2], where there are also three regimes, and the Stokes regime does not extend to r min — > 0. 

Having argued for the distinct identities of the inertially-limited-viscous regime versus the 
Stokes regime on theoretical grounds, we now offer evidence from experiment and simulation 



3 




min 



Figure 2: Phase diagram for 3D coalescence. (A) While the inertial regime [Tg 123 H31 13U IH01 HJ 122 SI EH1 EJ 
[SJ [55] and the Stokes regime [HJ [501 HH (HI 13U have been established in recent years, the inertially-limited- 
viscous regime is identified by this work. The inertially-limited-viscous to inertial crossover was recently 
determined by experiments (open circles) and a new scaling argument (solid line) |26j . Here, we identify 
the inertially-limited-viscous to Stokes crossover with a force-balance argument (dashed line: eqn. [1] with a 
proportionality constant of 1.4), simulations (filled triangles) and experiments (open triangles). The data 
departs from the prediction at large r m i„, where we expect finite size effects to enter. (B) To observe the 
inertially-limited-viscous to Stokes crossover, we measure Vb.o.d. versus r m i n /A over a range of viscosities from 
simulation (solid lines) and experiment (symbols). The velocities fall onto the creeping-flow curve at large 
rmin/A, but peel away at smaller r m i n /A. The velocity scaling at small r m i n /A is qualitatively captured 
by eqn. [2] The crossover neck radius at which the macroscopic drop velocity Vb.o.d. merges onto the Stokes 
solution is plotted in (A), for the viscosities that exhibit such a crossover within our range of data. 
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that these regimes are, in fact, different. First, we probe the global motion of the drops by 
measuring the velocity of the back of one drop, Vb.o.d.- Since this point is the farthest from 
the singularity, it isolates the global motion from the flow near the growing neck. 

For 3D drops in the inertially-limited- viscous regime, force-balance gives a c , , m . ~ F y /m = 
3 , yr m i n /2A 3 p. Using r min = ry/fj, as seen in Fig. [3j5 consistent with previous experiments 
E EES S 123 E8], we integrate to get: 

If Oh > 1, then the flows eventually enter the Stokes regime, where to first order: 
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We find that the creeping-flow simulation follows the exact analytic 2D Stokes solution 
(eqn. [3]). In Fig. [2f3, we plot Vb.o.d. for drops of finite viscosity in the simulation and exper- 
iment. The curves exhibit the predicted super-linear growth of Vb.o.d. a t early times until 
the velocities merge onto the Stokes curve. The data show exceptional agreement between 
simulation and experiment. 

We measure the crossover neck radius at which the macroscopic drop velocity Vb.o.d. merges 
onto the Stokes solution. We plot Oh versus r m i n /A as the threshold for entering the Stokes 
regime on the phase diagram for 3D coalescence (Fig. |2]A.). At higher viscosities, the linear 
(Stokes) regime is entered at smaller r min . The data agree well with the prediction from eqn. 

m 

The neck shapes also differ between the Stokes and the intertially-limited- viscous regimes. 
Because the exact analytic 2D solution extends over the entire domain of r m i n (whereas the 
3D solution only exists for small r min ), we use it here as a convenient way to account 
approximately for finite-size effects and to compare our experiment and simulation. In Fig. 
|3]A. and B, we plot this 2D Stokes solution [TJ] in the neck region against experiment and 
simulation, for Oh ^> 1 and Oh = 0.6, in the r-z plane (with the origin at the initial point 
of contact). The Stokes solution agrees with the high- viscosity data, but it clearly fails to fit 
the shapes at Oh = 0.6 where both experiment and simulation show a much broader neck. 
In particular, the 2D analytic Stokes solution has a maximum neck curvature, k, that obeys 
1/kA = \{r m i n j A) 3 to first order, which we find is in good agreement with our data in the 
Stokes regime (i.e. Oh 3> 1). In the inertially-limited- viscous regime (i.e. Oh ~ 1), we also 
find that 1/kA scales as (r min /A) 3 , but with a significantly larger prefactor (approximately 
1.2 instead of ~). 

The dynamics also differ between these two regimes. For high-viscosity drops, measure- 
ments of the neck radius r min versus time in the experiment and simulation are consistent 
with the exact analytic 2D Stokes solution (Fig. EP)Jll For lower- viscosity drops, the 2D 



1 Previous viscous coalescence experiments [32j EQl [TJ [33l [26] have compared the neck radius r m j n (r) 
against the theoretical prediction |15j . r m i n sa tj\ \ti(tj/ [iA)\/Ttfi. This approximate form breaks down |15j 
for r m i n > 0.03 A. Therefore, we compare our measurements against the full analytic 2D solution |19] . which 
can be done over the entire range of data. 
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Figure 3: Coalescence dynamics in the Stokes regime (left column: creeping flow simulation, Oh = 440 
experiment) versus the inertially-limitcd-viscous regime (middle column: Oh = 0.6 simulation and experi- 
ment) and the inertial regime (right column: Oh = 0.007 simulation and experiment). (A,B,C) Neck profiles 
from simulation (solid lines) and experiment (symbols) at three different times, compared with the 2D Stokes 
theory in (A) and (B) (dashed lines [H]). The Stokes profiles agree with high- viscosity data (A) but do 
not capture the broader interfacial shapes at Oh = 0.6 (B). (D,E,F) r m i n /A versus rescaled time in the 
simulation (solid lines) and experiment (circles). High- viscosity drops (D) follow with the Stokes theory 
(dashed line). At intermediate viscosity (E), r min /A does not agree with the Stokes theory but instead 
grows at the viscous-capillary velocity, r m i n = Tj/fj, (dotted line). At low viscosity in the inertial regime 
(F), r min /A = 1.4r 1 / 2 ( 7 / / o^ 3 ) 1 / 4 (dash-dot line). Note that in (E), we begin plotting the simulation data 
at tj//j,A — 6 x 10~ 3 , where we estimate that transients from the initial conditions have decayed. The 
experimental data in (E) was obtained by a high-speed electrical method on glycerol-salt-water drops (p, = 
230 cP, 7 = 65 dyn/cm, p = 1.2 g/cm 3 , and A = 0.2 cm) [26]. (G,H,I) Instantaneous streamlines from 
simulations at r m i n = 0.03A in the (G) Stokes regime (creeping flow simulation), (H) the inertially-limitcd- 
viscous regime (Oh = 0.6), and (I) the inertial regime (Oh = 0.007). The flows are qualitatively different in 
all three regimes. 
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Stokes solution does not fit the data (Fig. [3)5) . Instead, the neck radius grows linearly with 
time, consistent with r min = Tj/fj,, a form that one might guess from dimensional analysis 
alone. (For the experimental data in Fig. [313, we coalesce hemispherical drops attached to 
circular nozzles separated by a distance 2A. This altered boundary condition does not affect 
our results: using high-speed imaging, we find that for r m j„ <C A, the dynamics is insensi- 
tive to this change in boundary conditions in both the inertially-limited-viscous and Stokes 
regimes.) This linear growth has been observed in previous experiments, but has incorrectly 
been assumed to be the dynamics of Stokes coalescence [U EDI SI ESI E3] , and therefore was 
not recognized as evidence of a new regime. We emphasize that the observed power-law, 
Tmin oc t, is different from the inertial scaling where r m j n oc t^' ^ (as shown in Fig. [3]F and by 
previous work [Ii2[25j[l3j[3T]), which demonstrates that the inertially-limited-viscous regime 
is distinct from the inertial regime as well. 

Lastly, our simulations give the flow profiles near the singularity in the Stokes and the 
inertially-limited-viscous regimes (Fig. [3p and H). The flow is expected to occur over a 
length scale comparable to r min in the Stokes regime [15] . and over a length scale r^ in /A in 
the inertially-limited-viscous regime [26]. Indeed, while features in the creeping flow stream- 
lines are roughly the size of r min , the streamlines at intermediate Oh exhibit recirculation 
zones, which constrict the flows near the neck. Comparing the geometry of the streamlines 
further solidifies that the new regime (Fig. (3j3,E,H) is distinct from the one described by 
pure Stokes flow (Fig. |3A,D,G) and the one described by inertial flow (Fig. |3C,F,I). The 
streamlines in Fig. [3J corroborate the definition of the Reynolds number proposed in [26] 
that dictates the inertially-limited-viscous to inertial crossover. 

Conclusion 

As two drops begin to coalesce and a microscopic liquid neck forms between them, the 
curvature of the interface and the Laplace pressure that develops due to surface tension 
both diverge at the instant when the drops first touch. In drop coalescence with no external 
fluid, previous work [T5] incorrectly led to the conclusion that only viscous forces, along 
with surface tension, should dominate on small scales, a dynamical regime referred to as 
the Stokes regime. Our work identifies a necessary condition for Stokes flow to occur. The 
dynamics cannot be in the Stokes regime until the surface tension force around the neck is 
large enough to rigidly translate the two initially stationary drops towards each other. The 
inexorable resistance of inertia rears its head at even these small scales. 

Therefore, a heretofore unknown dynamical regime controls the singularity at early times 
for drops of any viscosity. In this initial, asymptotic regime of drop coalescence in air, all 
of the underlying forces, that is inertial, viscous, and surface tension forces, are important. 
Hence, the two dynamical regimes referred to as the Stokes regime, where inertia is negligible, 
and the inertial regime, where viscous force is negligible, can only be attained once the neck 
has grown to a sufficient size. Once the drop has entered the Stokes or inertial regimes, our 
measurements are consistent with the earlier predictions for the dynamics in those regimes. 

A dynamically similar response is observed when a liquid filament breaks in air. At small 
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neck size, the viscous and inertial regimes both give way to a third regime, where inertial, 
viscous and surface tension forces are all important [2]. When a liquid filament breaks in 
another liquid, however, the dynamics are qualitatively different: in that case the asymp- 
totic dynamics of thinning may occur in the absence of inertia [231 E]- Further insight may 
likewise be gained by studying drop coalescence inside a second immiscible liquid, which is 
a problem of immense practical importance [HI EE] • 

Materials and Methods 

The experiments were conducted by authors JDP, JCB and SRN. High-speed imaging and 
electrical measurements were separately performed to capture the coalescence dynamics of 
isolated liquid drops in air. In high-speed imaging measurements, two silicone-oil drops with 
radii i ^ 0.1 cm were suspended side by side as in Fig. [TA from syringe needles. Different 
silicone oils were used in order to vary the liquid viscosity, /i, from 49 cP to 58,000 cP, while 
keeping other fluid parameters constant (surface tension, 7 = 20 dyn/cm and density, p = 
0.96 g/cm 3 to 0.98 g/cm 3 ). Because they are highly wetting, pendant silicone-oil drops tend 
to climb up stainless steel syringe needles until the needles protrude from the bottoms of the 
drops. To prevent this, the needles were treated with an electronic coating (Novec™ EGC- 
1700, 3M) that inhibits wetting. The drops were aligned and then slowly translated with 
a micrometer stage until they gently touched at their equators. The resulting coalescence 
dynamics were recorded with a high-speed digital camera (Phantom vl2, Vision Research). 

The electrical method is described in detail in [26]. The experimental data in Fig. 05 
was obtained by this method on glycerol-salt-water drops (fi = 230 cP, 7 = 65 dyn/cm, p = 
1.2 g/cm 3 , and A = 0.2 cm). 

The simulations were conducted by authors SA, MTH and OAB. The coalescence of two 
identical, isolated spherical drops of radius A of an incompressible Newtonian fluid that 
are surrounded by a dynamically passive gas is simulated by connecting them with a small 
bridge of radius r min (typically equal to 0.1% of the drop radius) and height z min <C r min . 
The ensuing coalescence dynamics are governed by the continuity and the Navier-Stokes 
equations, i.e. the Navier-Stokes system. Because the two drops are identical and the two- 
drop configuration is axially symmetric, the computational domain is the planar quadrant 
that consists of one of the drops and one half of the bridge that is bounded by the plane 
of symmetry, the axis of symmetry, and the liquid-gas (L-G) interface. The Navier-Stokes 
system is solved subject to symmetry boundary conditions along the plane of symmetry and 
the axis of symmetry, and the kinematic and traction boundary conditions along the L-G 
interface [TJ, [29] . This free boundary problem is solved numerically by a fully implicit method 
of lines (MOL) arbitrary Lagrangian-Eulerian (ALE) algorithm that uses the Galerkin/finite 
element method (G/FEM) for spatial discretization and an adaptive finite difference method 
(FDM) for time integration [71 |2H] . On account of the free boundary nature of the problem, 
the interior of the flow domain is discretized by an adaptive elliptic mesh generation algorithm 

The G/FEM converts the transient system of nonlinear partial differential equations 
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(PDEs) to a system of nonlinear ordinary differential equations (ODEs). The FDM time 
integrator reduces the system of ODEs to a large system of nonlinear algebraic equations. 
This system of equations is then solved by Newton's method with an analytically calculated 
Jacobian. 

Starting from an initially quiescent state, the dynamics are followed until the two drops 
have coalesced into one and the dynamics have ceased. Simulations are carried out for both 
situations in which inertia is present, i.e. Oh is finite, and also when inertia is negligible, i.e 
1/Oh = 0, such that the drops undergo creeping (Stokes) flow. 
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